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Abstract 



Description of the magnetic Bianchi VIo cosmologies of LeBlanc, Kerr, and 
Wainwright in the formalisms both of Belinskii, Khalatnikov, and Lifshitz, and 
of Misner allows qualitative understanding of the Mixmaster-like singularity 
in those models. 
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LeBlanc, Kerr, and Wainwright (LKW) have recently used the dynamical systems 
approach popularized by Wainwright and Hsu [fj] to study the class of Bianchi VI cosmolo- 
gies containing a magnetic field and a perfect fluid. As part of their complete analysis of 
this system, they reported that the generic singularity (for most fluid equations of state) 
was of the Mixmaster type [[5],|J. The purpose of this comment is to elucidate the nature 
of this feature of the models by reexpressing the results in the Belinskii, Khalatnikov, and 
Lifshitz (BKL) and Misner's minisuperspace (MSS) terminologies. While nothing new is 
thereby revealed, those more conversant with these formalisms will find LKWs somewhat 
surprising result less mysterious. In most of the following, we shall restrict attention to the 
electrovac case, only remarking at the end on the effect of inclusion of a perfect fluid. 

First consider the "standard Mixmaster model" (SMM) of (diagonal) vacuum, Bianchi 
IX. The (BKL) metric components a 2 , b 2 , and c 2 are related to the MSS variables Q and f3± 
describing volume and anisotropy by 

a = In a = Q — 2j3 + , 

( = \nb = + /3 + + V3/3-, 

7 = lnc = n + (3 + -V3(3- (1) 

Einstein's equations may be obtained by variation of the superhamiltonian A/"7i[] with 

2H = -p 2 n +p 2 + +p 2 _ + U(n,f3±) (2) 

where p^ and p± are respectively canonically conjugate to Q and (3± and 

U = e m V((3 +} f3.) 

= e 4« + e 4C + e 4 7 _ 2e 2 (^) - 2e 2(a+7) - 2e 2(a+c) . (3) 

Arbitrary constants in front of U may be absorbed by a suitable choice of spatial coor- 
dinates. In (|J), the MSS potential, U, is related to the spatial scalar curvature through 

x The rescaled lapse is = Afe~ 3Q where dt = Mdr with r our choice of time coordinate and t 
comoving proper time. In LKW, • = d/dt while M = 1 yields the BKL time coordinate. 
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U = e 6n 3 R where e 6n is the determinant of the spatial metric with scalar curvature 3 R. 
Direct comparison of (fj) and (Q) with the LKW form of the Hamiltonian constraint suggests 
the identification 

2a 

S± = ^ ; N 1 = c l — , (4) 
-pn -pn 

etc. LKW's equations of motion are then obtained with J\f = — p^ 1 so that the time 
coordinate is fl The c% contain factors from the structure constants. Thus change of 
Bianchi type can lead to vanishing of some of the iVj and sign changes of others @. 

The approach to the singularity (Q — ► — oo) in the SMM has long been known to be 
characterized by an infinite sequence of Kasner models (vacuum Bianchi I) described by the 
single BKL parameter u @,f|]. The Kasner solution is just that obtained for U = in ([|). 
The superhamiltonian becomes that for a free particle so that the momenta pq and p± are 
constant. The Kasner epoch changes when a scattering (bounce) off the MSS potential U 
occurs.0 The values of u in successive Kasner eras are related by the BKL map ||, equivalent 
to the Kasner map described by LKW [|TJ. All discussions of chaos in the SMM derive from 
the well-known sensitivity to initial conditions of this map (see for example ||). Here we 
need consider only the features of the SMM which are required to derive the BKL (or Kasner) 
map since the map will result in any model in which these same properties are present. In 
fact, only the first three terms in U are needed.^ Since the Kasner solution is characterized 
(in the approach toward the singularity) by two contracting metric components and one 
expanding metric component, one of a, (, or 7 will dominate. For convenience, choose the 
dominant one to be a so that U ~ e 4n-8/3+_ canon i ca ] transformation generated by 



2 The overall e 4 ^ dependence of U means that, as ft — > —cc, U — > 0. However, the Kasner solution 
(expressible as f3± = 0± + E±(fi — fio)) allows terms in U to be of order unity (e.g. if Q = 2/3+). 

3 One can argue that the remaining terms are always dominated by one of the first three except 
near the 120° corners of the MSS potential In this discussion of generic Mixmaster behavior, 
the "anomalous" corner behavior will be ignored. See, however, |3]||. 



x = Q — 2(3+ and y = —2Q + f3 + puts the Hamiltonian constraint @ into the form 

Spl-3pl+p 2 _ + e^ = 0. (5) 

Eq. (|5|) describes scattering off a potential with relationships among the momenta before 
and after the scattering of constant p y and p_ with p x — > —p x . Since the BKL parameter 
u may be expressed in terms of these momenta ||, the scattering rules lead immediately to 
the BKL map. It is important to note that, even though the exact solution corresponding 
to (H) is known, only the momentum rules are used. 

With the identifications (§]), the Hamiltonian constraint yields the identity 

p^U = 1 - S (6) 

(where £ = E 2 + + E 2 . 0) so that 

^lnp, = 2-g (7) 

where q = 2S. The LKW equations for the N{ then yield the relationship of these variables 
to the MSS or BKL ones. For example, N% = rii/pQ in |T|] 

^=fa-™ (8) 

with (0) has the solution (using d(3±/dVL = £±) 

n 1= n?e 2 ^ + . (9) 

Similar constructions can be performed for N2 and N3. 

The most attractive feature of the LKW formalism is its ability to handle all the Bianchi 
types with the same equations. Thus extension to the magnetic (but otherwise empty) 
Bianchi VIo equations is straightforward. To make the connection to the BKL and MSS 
formulations, however, we must include some properties of the structure constants (i.e. in 
the Cj of (|P). Thus, for Bianchi VIo, c± — and C2C3 < while Hi = hi/pn is the only 
magnetic field component consistent with this model ffl. The magnetic field contribution 



to the Hamiltonian constraint leads to a modification of (|6]) such that p^ 2 U = 1 — £ — |i/f 
and q = 2S + f-^i • These identifications allow the LKW equation for Hi JIJ 



(q-l-2X+)H 1 (10) 



to yield the solution 



h x = h\e a - 2 ^. (11) 
It is now possible to transcribe the LKW constraint for this model into the MSS language 



as 



m = -pi +p 2 + + pi + [e^ + e 4 ^ + 2e 2 ^] + £e 2a . (12) 

We note that, compared to the SMM constraint (0), the Kasner (momentum) part is un- 
changed. The first two terms in the brackets, identical to those in SMM, are two of those 
required for the BKL map derivation. The third term is the irrelevant cross term (with 
irrelevant changed sign). The final term in (|T2"D, from the magnetic field, is essentially the 
square root of the corresponding term in the SMM potential. However, the arbitrary con- 
stants, the absent cross terms, and the difference in the power of e a have no effect on the 
BKL map derivation. The same canonical transformation to x and y can be made yielding 
a model for the scattering from one Kasner epoch to another of 

3^-3^+p^ + £e 2 * = (13) 

rather than (|5|). Even though (as pointed out by LKW) the exact solution has changed, the 
rules describing the change in momenta at the bounce have not.[] Thus the BKL map may 
be derived as a property of the magnetic Bianchi VIo vacuum model. 

Finally, we note that a perfect fluid will enter the Hamiltonian constraint (0) as a term 
U — > U + poe 3 ( 2 ~ r ^ n for equation of state p = (T — l)p. As — > — oo, this term goes to zero 



l Of course, (||) is exactly reproduced to describe bounces for either £ or 7 approximately zero. 
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(if T 7^ 2 in the range discussed by LKW) and can never revive due to cancellations such 
as those in other terms of U that produce Mixmaster bounces. Only if T = 2 will the fluid 
play a role in the singularity approach (as has been found by LKW). 

Thus, the connection between the LKW formalism and the BKL and MSS descriptions 
of Bianchi cosmologies clarifies for those familiar with the latter the similarity between the 
singularity dynamics of the SMM and magnetic Bianchi VIo models by showing that the 
differences between the two fail to invalidate the derivation of the BKL map. Of course, the 
BKL approximate description of the SMM contains more than just the parameter u ||[7|]. 
While the exact solution to the bounce model equation is still not needed for this expanded 
description, differences between the two models might prove interesting. 
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